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The development of one-dimensional parametric instabilities of intense long-wave 
plasma waves is considered in terms of the so-called hybrid models, when electrons 
are treated as a fluid and ions are regarded as particles. The analysis is performed for 
both cases when the average plasma held energy is lower (Zakharov’s hybrid model - 
ZHM) or greater (Shin’s hybrid model - SHM) than the plasma thermal energy. The 
efficiency of energy transfer to ions and to ion perturbations under the development 
of the instability is considered for various values of electron-to-ion mass ratios. The 
energy of low-frequency (LF) oscillations (ion-sound waves) is found to be much 
lower than the hnal ion kinetic energy. We also discuss the inhuence of the changes 
in the damping rate of the high-frequency (HF) held on the instability development. 
Reduced absorption of the HF held leads to the retardation of the HF held burnout 
within plasma density cavities and to the broadening of the HF spectrum. At the 
same time, the ion velocity distribution tends to the normal distribution in both 
ZHM and SHM. 
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I. INTRODUCTION 


Powerful fields that can be excited and are able to propagate in a plasma, i.e. in a 
nonlinear medium with a wide range of characteristic oscillations, often produce a great 
variety of instabilities. As a result, the new generally unsteady states come into existence in 
the medium. The theoretical description of all these successive relaxation processes in the 
form of cascades of instabilities requires solving very complex problems (in the mathematical 
sense) with the use of numerical methods and numerical experiments. 

Investigating the after-effect of instabilities is important for a variety of applications 
in plasma physics and plasma electronics. Moreover, the employment of more and more 
powerful sources of energy leads in particular to the use of plasma in the conventional 
vacuum high-current electronic devices. The experimental and practical activities in space 
must also take into account the presence of plasmas. The very great variety of instabilities 
and even cascades of instabilities create problems for the implementation of the technologies 
developed up to now. On the other hand, the abundance of various unstable states provides 
opportunities for discovering and developing entirely new technologies. 

Among the variety of plasma instabilities, the class of parametric instabilities should be 
especially emphasized. The modulation of velocities and densities of plasma particles by 
intense external helds results in the development of instabilities similar to those described 
by Mathieu’s or Hill’s equations. In the case when the external radiation is in resonance 
with plasma eigenmodes, such instabilities can be interpreted as decay or modulation ones. 
The interest to the parametric instability of intense Langmuir waves, which can be easily 
excited in the plasma by various sources^”- was stipulated, in particular, by new possibilities 
in heating electrons and ions in a plasma. The correct methods for the description of the 
parametric instability of long-wave plasma waves were developed in the pioneer works of 
V.P. Silirt^ and V.E. Zakharov^^^— . 

The theoretical concepts proposed by V.P. Silin^- were conhrmed by the early numer¬ 
ical experiments on the one-dimensional simulation of the parametric decay of plasma 
oscillations!^ (see also article!^'!^, and review!^). However, the greatest experimenters’ inter¬ 
est was provoked by the mechanism of wave-energy dissipation discovered by V.E. Zakharov. 
The analytical studies, laboratory-based experiments, and numerical simulations, performed 
at an early stage of studying these phenomenap^^— , have conhrmed the fact that in some 
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cases, in the course of instability, a significant part of the pump-held energy is transformed 
into the energy of short-wave Langmuir oscillations2i>2^ accompanied with bursts of fast 
particlesl^”— . The modulation instability of intense Langmuir waves in non-isothermal plas¬ 
mas also leads to collective ion perturbations, in particular, to the generation of ion-sound 
waves2^^— . 

In Zakharov’s model^ that describes the instability of intense long-wave Langmuir waves 
in a non-isothermal plasma, namely the modulation instability, results in the excitation of a 
range of short-wave oscillations. In Shin’s model, a strong Langmuir wave in a cold plasma 
leads to intense oscillations of the electron velocity with the amplitudes comparable to the 
wavelengths of the excited modes. In this case, the instability in the general case should 
be referred as parametrio^. Nevertheless, both these processes are similar— Hence, the 
term “modulation instability” can be applied to the instability of a strong Langmuir held 
within the framework of Silin’s model. 

In the works^i^ an attempt was made to compare these models, which have similar 
physical nature, by the example of one-dimensional description. The choice of the one¬ 
dimensional approach, as was noted by J. Dawson^, “often keeps the main features of the 
processes, but simplihes their description and leads to a fuller understanding of what the 
important phenomena are” (see also^). Of particular interest is the process of ion heating, 
so in this paper we use the particle (or hnite-sized particle) description for ions because the 
account of inertial effects can be signihcant just at the nonlinear stage of the process^. It 
was observed^i^ that simulation in terms of the so-called hybrid model (incorporating one of 
Zakharov’s equations for the HF waves and using particle simulation for ions) demonstrates 
that fluctuations of ion density are rather signihcant and favor the development of parametric 
instability. The non-resonant interaction between ions and HF plasma oscillations, along 
with ion trapping by the potential wells produced by these oscillations, lead to an instability 
of density cavities resulting from the modulation instability and produce fast particle groups. 

Clark et ah— compared the hybrid model with Zakharov’s hydrodynamic model. Due to 
higher level of ion density huctuations the number of cavities in the hybrid model is found 
to be signihcantly greater than in Zakharov’s model and their depth is smaller—. Inte¬ 
gral characteristics of the instability after-effects for both models turn out to be essentially 
identical. 

In papers^^i^ the simulation of one-dimensional ion dynamics was performed in terms of 
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the particle method^>^. The number of particles used in numerical calculations was 2 • 10*^, 
which is equivalent to the number of ions, about (2 ■ 10^)^ ~ 10^^, in the three-dimensional 
case, in agreement with the conditions of most experiments. Thus, the interaction between 
modeling particles and plasma oscillations in this simulation is in rather good accordance 
with the interaction between real particles and plasma waves, naturally with regard to the 
inherent limitations of the one-dimensional description. Nevertheless, there is reason to 
believe that the description of held energy transfer to ions within the framework of the 
hybrid model represents the real conditions of ion heating by intense Langmuir oscillations 
in plasmas. Moreover, the one-dimensional description makes it possible to select arbitrary 
electron-to-ion mass ratios. 

Below, we discuss the efficiency of energy transfer from Langmuir oscillations to ions and 
ion density perturbations under the development of the modulation instability in both cases 
of non-isothermic hot and cold one-dimensional plasmas within the framework of hybrid 
models and for different values of the electron-to-ion mass ratio. The attention is mainly 
concentrated on the effect of HF held burnout within density cavities accompanied by energy 
transfer to the ion component of the plasma. 

II. THE HYBRID MODELS (HM) OF PARAMETRIC INSTABILITY 
A. The hybrid model based on Shin’s equations (SHM) 

When the intensity of the external electric held is much greater than the specihc thermal 
energy of plasma electrons W = \Eq\^/Att 3> noTg, it is reasonable to explore the approach 
presented by V.P. Silin^. 

Let consider a one-dimensional plasma system where an intense plasma wave with the 
wavelength Aq and frequency Uq is excited by an external source. This intense wave will be 
referred to as the pumping wave. Since the parametric instability results in the growth of 
oscillations with rather small wavelength A -C Aq, the pumping wave can be suggested to be 
spatially uniform within the region of interaction: 

Eq = —i{\Eo\ exp{ia;ot + i4>} — |Eo| exp{—iwof — i(j)})/2, (1) 

where |Eo|and 0 are the slowly varying wave amplitude and the phase respectively, cjq is 
the external wave frequency, Uq and Te are the density and temperature of plasma electrons. 
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Charged particles of plasma oscillate under the action of the electric held and their velocities 
can be written as uoa = — (e^ |-Eo| coscf) = —cjo&o cos0, where bo = Cq, |i?o| /fnaUJ^ is 

the particle oscillation amplitude. 

The equations, governing the nonlinear dynamics of the parametric instability of an 
intense plasma wave, were derived in the work—. The equations for the HF plasma held 
spectrum modes E = ■ exp{inkox) (plasma electrons are treated as huid and 

described by hydrodynamic equations) are given by 


dEr, 

dt 


cjpg — (Ug n® 

- -En + e — E„ 

2uo nlj 
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— '^ 2^ — 0. (2) 


Here Upe = is the background electron plasma frequency, e and rUg are the 

electron mass and charge, M is the ion mass of an ion, En = |i?n| •exp(i'^„) is a slowly varying 
complex amplitude of the electric held of electron plasma oscillations whose wave-number 
is kn = nko, ko = 2n/L, where L is the characteristic dimension of the plasma system, 
~ ' exp(m/co^) IS the ion charge density, an = a ■ n, n,m are nonzero integers, 

and ±1, e.g. a„ = nkob = n (e/coF^o/'^e<^o)- The term En6n^/n\j in Eq. ([2]) simulates the 
HF plasma wave damping by electrons, with hm = 20. 

The equations of ion motion are given by 

^ exp{i/cona;4, (3) 


and the ion density can be determined from 
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The slowly varying electric held strength En, acting on the ions, is equal to 


En = 


Att? 2 
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The equation for the uniform component of the electric field Eq = |ii^o|exp(i0) can be 
written as 

pijp 

Y. ^r,-m[E*_^J2{am)e^^^ + E^Jo{aJ]. ( 6 ) 

^ m 

Note, that the values with subscripts with unlike signs are independent. In Eqs. ©-([HD, 
we have used the formulai^ 

OO 

exp{zasin$} = Jm{a) exp{im^}, (7) 

m=—OO 


where Jm{x) is the Bessel function. 

The normalized frequency shift A = {Up^—ojQ)/2Supe reaches the value of (me/2M)^/^ 
for the mode with the maximum growth rate of the parametric instability 6 (see^) 


5/cjpe 
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B. The hybrid model based on Zakharov’s equations (ZHM) 


As is shown in^S, Eqs. ©-© reproduce the equations obtained in^I after following sub¬ 
stitutions: (cUpg — a;Q)/2a;o (ai^g — uJq + kln?v^^)/2 ujq and Eq —iE^, E^ —)■ iE^ under 
the condition a„ 1, that means that Ji(a„) ~ an/2, Jo{an) ~ 1, ^ 2 (an) ~ ®n/8 


dt 


2a;o 2no 


m^O 


= 0 . 


(9) 


Note that the damping of the HE modes (the term oc 9) was introduced in Eq. © phe¬ 
nomenologically, because we are only interested in how the changes in the absorption level 
influence the ion heating. More precise analysis of the effect of kinetic processes on the 
Landau damping and on the saturation of parametric instabilities in an electromagnetically 
driven plasma can be found in papers^*^. 

In this case, the slowly varying electric field amplitude takes the form 


En — 


ikone 

AmUp 


f EnEQ + EqE*_^ + 'Y En-mElm j , 

\ m^O^n / 


that makes it possible to describe ions as particles within the context of of Eqs. 
pump wave amplitude Eq is governed by the equation 

OEq . Uq 


dt 2n( 


^ ^ 'ai^—mEm 0. 


( 10 ) 
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In this case, the growth rate of the parametric instability normalized to the plasma 
frequency is given by2^ 



( 12 ) 


III. STATEMENT OF THE PROBLEM AND THE INITIAL CONDITIONS 

The purpose of this paper is to clarify the efficiency of the energy transfer to ions and ion 
density perturbations in the course of development of modulation instability for the cases of 
both non-isothermal hot and cold plasmas in terms of the hybrid models. Both SHM and 
ZHM were considered for two cases of light and heavy ions.The parameters of the simulation 
are presented in Table 1. It is also interesting to elucidate the effect of HF spectrum damping 
and subsequent burnout of the Langmuir field within density cavities on the energy transfer 
to plasma ions. 

Below we employ, unless otherwise specified in the text, the following initial conditions 
and parameters. The number of particles simulating the dynamics of ions is 0 < s < S = 
20000. The particles are distributed uniformly over the interval —l/2<^< 1/2, ^ = 
kox/271, initial ion velocities are defined as d^s/dT\r=o = t’s|r=o = 0, the number of spectrum 
modes is —N < n < N, N = F/IOO. The initial normalized amplitude of the pumping 
wave is ao(0) = ekoEQ^O)/meOJp^ = 0.06,. The initial amplitudes of HF plasma oscillations 
are defined by the expression en\r=o = ^no = (2 + ffn) ■ 10“^ for the Silin model and by the 


TABLE 1; Simulation parameters for the hybrid models 


Model 


Light ions 


Heavy ions 
me/M = 8 • 10“® 


M/nie = 2 • 10^ 


SHM 


(me/M)(w^/(5^) = 0.43 
6/uJo = 0.44 • = 0.034 

UIq/S ~ U!pel5 = 29.4 


{me/M)[ujl/5'^) = 0.1 
5/uJo = 0.44 • = 0.0088 

Wq/iJ ~ Wpe/lJ = 113.6 


ZHM 


= 2noElW = 20 

ujq/5 = 2 = 282.6 

^/wq = SjlOpe = 3.5 • 10“^ 


{me/M){ujl/6^) = 2noTjW = 20 
wo/J = 2(noTe/lF)F2(M/TOe)F2 = 2234.4 
5/ujq = Sjujpe = 4.5 • 10“'^ 
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expression en\T=o = ^no = (0.5 + gn) ■ 10“"^ for the Zakharov model, where gn G [0; 1] is a 
random value, ekoEn/meUp^ = Cn exp^i'ijjn)■ The initial phases of spectral modes '0n|T=o are 
also randomly distributed in the interval 0-G27r. The ion density fluctuations and slowly 
varying electric held En are described by the dimensionless representations 

/ niko 

exp (2 

■n jko 

and 

ck^E^!TYl^UJp^ ^nr T 

The development of the instability was considered in terms of hybrid models Eqs. ()2|)- 
([H]) and Eqs. in our previous papers^^i^E. Here we give some results. The rate of 

damping of HE modes governs the rate of the held energy burnout in density caverns, from 
where the HE held has forced out charged particles The main part of the instability energy is 
initially concentrated in the HE Langmuir oscillations in parallel with the formation of the LF 
spectrum of density perturbations. Then the energy of the HE spectrum is transferred mainly 
to electrons. Thus, the shaped density cavities collapse, the trajectories of ions intercross, 
ion density perturbations become smoother and their characteristic scale growths with time. 
The relationship between ionic perturbations and the HE held is weakened and the instability 
is saturated. The amplitude of the main wave stabilizes after several oscillations at rather 
low level. The bulk energy is now contained in the perturbations of the electron component 
of the plasma. Some small portion of the initial energy transforms into the kinetic energy 
of ions. The estimate of the energy density transmitted to ions E^i^ can be obtained from 
the expression 


E, 


kin 


ITn 


0.27./.L.il, 

• 


m oj, 


(13) 


pe 


where Wq is the initial energy density of the intense Langmuir wave, I = Y2s i.d^s/d'T'Y 
normalized ion kinetic energy and 5 is the rate of the linear instability. The portion of energy 
transferred from the intense Langmuir wave to ions is determined by the ratio hho/'^oT'e for 
the case of non-isothermic plasma (ZHM) and by the ratio for the case of cold 

plasma (SHM). 

Below we consider more closely the nature of the energy redistribution with time and 
especially the process of energy transfer to the LF perturbations. We also discuss the 





(a) ZHM, light ions (b) SHM, light ions 




(c) ZHM, heavy ions (d) SHM, heavy ions 

FIG. 1; Time evolution of relative values of: energy of the main Langmuir wave (1), energy 
of low-scale plasma wave spectrum (2), energy transferred to electrons (3) and ions (4). 

specific features of the excitation of LF ion-sound waves in both non-isothermic and cold 
plasmas. More attention will be focused on the role of absorption of HF spectrum that is 
responsible for the burnout of the HF held in the density cavities. We investigate the effect 
of this process on the excitation of the LF spectrum and most importantly on the kind of 
ion velocity distribution function and on the proportion of the total energy transferred to 
ions. 

IV. THE RESULTS OF NUMERICAL SIMULATION 

Figured] shows the energy redistribution between the main Langmuir wave, the small-scale 
plasma wave spectrum and plasma electrons and ions for the following values of parameters 
determining the damping rate of the HF spectrum: hm = 20, 0 = 0.005. 

The analysis of the numerical simulation results shows that the energy of intense long¬ 
wave Langmuir waves is hrst transferred to short-wave Langmuir oscillations. Just at this 
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(c) ZHM, heavy ions (d) SHM, heavy ions 

FIG. 2: Dependence of the amplitude of the LF modes Mn (1) and the frequency d^n/dr 
(2) on the wave-number at the stage of developed instability. 

stage the cavities of plasma density, hlled with HF plasma oscillations, are formed. After 
that, the HF held burns out due to the damping on electrons that is included in the hybrid 
models phenomenologically. The energy of the HF held therewith converts into the energy 
of plasma electrons. Under these conditions, the cavities collapse and thus excite LF waves, 
the ion trajectories intercross, and the energy of both collapsed caverns and LF spectrum is 
transferred to ions. 

The root-mean-square velocity of ions, (7{v) = at the hnal stage of the nu¬ 

merical simulation is equal to a{v) = 0.015 for the case of light ions and a{v) = 0.006 for 
heavy ions in ZHM and, respectively, to a{v) = 0.002 for light ions and a{v) = 0.0005 for 
heavy ions in SHM. The total kinetic energy of ions in assumed units I = Yls {^^s/drY is 
equal to 4.689 for the case of light ions and 0.808 for heavy ions in ZHM and a(v) = 0.086 
for light ions and 0.005 for heavy ions in SHM. The variations in the values of the total 
energy are caused by different linear growth rates in the two models under consideration, 
and by different ion masses in the simulation of light and heavy ions. The hnal ion velocity 
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(a) ZHM, light ions 



0.2 0.4 



(b) SHM, light ions 



(c) ZHM, heavy ions (d) SHM, heavy ions 

FIG. 3: Evolution of the ion kinetic energy (1) and LF field energy multiplied by factor 70 

(2) with time. 

distribution can be fitted by the normal curve with the use of the values of rms velocity. 
The particles outside the normal distribution (mainly in the so-called ’’tails”) possess 13.8% 
of the total energy for light ions and 9,2% for heavy ions in ZHM model and much more in 
SHM: 25,6% for heavy ions and 13% for light ions, respectively. It means that in the case 
of instability of the intense wave in a cold plasma, a signihcantly greater proportion of fast 
particles should be expected. 

We are interested not only in the ion kinetic energy distribution, but also in the collective 
excitation of ion oscillation (see Fig. |2]), hence we dehne the frequency of the mode with the 
wave vector nko associated with these oscillations, i.e., 

_ d_ f Mnr ^ j ( 

dT ^ dT ( + Ml, j / \ VMJ. + Aft 

where the phases of LF modes can be found from the expression 

Mn = Mnr + iMni = \J • • exp(i<I)„). 
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(e) ZHM, 0 = 0.001 


(f) SHM, 0 = 0.001 


FIG. 4: The ion velocity distribution for the case of light ions. 


It should be noticed that the intensity of the LF spectrum in the case of a non-isothermal 
plasma (ZHM) is quite high in a wide range of wave numbers, that corresponds to the 
spectrum of ion sound after the destruction of density cavities detected in numerical 
experiments^. In a cold plasma (SHM), in contrast, the long-wave oscillations dominate in 
the spectrum. 
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For both models, the ion kinetic energy in the assumed units can be written as 


1 

2 



- 1/2 



( 15 ) 


and the energies of collective excitations for ZHM and SHM, respectively, reduce to 


1 m 1 6 

Stt^ M ujpe 




1 m 








(16) 


Note that in Zakharov’s model these oscillations are referred to as ion-sound waves. 

Figure [3] demonstrates the time evolution of the ion kinetic energy and LF held energy. 
It should be noted that the energy of the LF held is far smaller than ion energies in all cases. 
Reducing the held energy with time is caused by the energy transfer to ions as well as by 
the destruction of plasma density cavities,—. 

The rate of the HF held burnout within density cavities is determined by the value 
Q = 6/6. It is of interest how the simulation results depend on this parameter. Obviously, 
the decrease of this parameter not only inhibits the burnout of the HF held in the cavities, 
but also broadens the spectrum of HF modes, i.e. it increases the contribution of small-scale 
components that leads to the deepening of plasma density cavities and to the growth of the 
kinetic energy of ions ejected from the cavities. 

Note that for both models the ion velocity distribution function approaches the Maxwellian 
curve with decreasing damping rate of HF modes, as may be seen from Fig. 01 

Table m demonstrates the extent of deviation of the ion velocity distribution function, 
obtained by numerical modeling, from the htted Maxwellian curve for the cases shown in 
Fig. dH). 


TABLE H: Deviation of the ion velocity distribution function, obtained by numerical 
simulation, from the htted Maxwellian curve 


Damping rate 

ZHM 

SHM 

0 = 0.05 

19.9% 

13% 

0 = 0.015 

9.9% 

13.4% 

0 = 0.001 

6.9% 

8.8% 
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(a) ZHM (b) SHM 

FIG. 5: Time evolution of the LF spectrum energy for the case of light ions, 
1 - 0 = 0.05, 2 - 0 = 0.015, 3 - 0 = 0.001. 




(a) ZHM (b) SHM 

FIG. 6: Time evolution of the ion kinetic energy for the case of light ions, 

1 - 0 = 0.05, 2 - 0 = 0.015, 3 - 0 = 0.001. 

Figure ([5]) shows that in the case of a non-isothermal plasma the maximum energy of ion- 
sound oscillations remains practically unchanged as the damping rate of HF held decreases, 
whereas the formation of the LF spectrum occurs with higher rate. In a cold plasma, on 
the contrary, the intensity of LF oscillations grows with the decrease of the damping rate 
of the HF held. After that, the LF spectrum is suppressed and its energy is transferred to 
ions. As might be expected, the energy, transmitted to ions, increases with the decrease of 
the damping rate of HF oscillations almost in the same proportion in both non-isothermal 
and cold plasmas (see. Fig. (j6])). 

It shonld be noted in conclusion that the ion-density pertnrbations with spatial scale 
smaller than the Debye radius rj^i = VTi/oOpi do not contribute to the formation of LF 
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electric fields by virtue of the screening effect. The Debye radius can be estimated from the 
expressioii^ 


rDiko/2n 


Rdi OC 


Vjko 

271^L 




1/2 


= < Vs > 




(17) 


At the stage of developed instability this value is of the order of Rm < 10“^ and the 
number of spatial spectral modes of ion density does not exceed I/Rdi that is in agreement 
with the previous analysis. 


V. CONCLUSION 

The most important consequence of the instability of intense Langmuir waves in plasmas 
is the transfer of a portion of the held energy to ions and LF plasma density oscillations. It is 
reasonable to consider this problem within the framework of hybrid models, where electrons 
are treated as a huid and ions are regarded as particles. The dynamics of the instability of 
intense long-wave oscillations in both hot and cold plasmas appears to be similar— >21. 

It was noted earlier— that in the case of hot plasmas ions acquire a portion of held energy 
that is proportional to the ratio of the held energy to the plasma thermal energy. In the 
case of cold plasmas, ions get a portion of held energy that is proportional to the ratio of 
the instability growth rate to the plasma frequency or that is the same as to the cubic root 
of the ratio of electron to ion masses. The energy transferred to ions in the case of heavy 
ions is signihcantly smaller than in the case of light ions. Moreover, the amount of energy 
transferred to ions in the case of the cold plasma is inversely proportional to the cube root 
of the ion mass, and in the case of the hot plasma the portion of energy transferred to ions 
decreases with growth of the ion mass faster— >21. The kinetic energy distribution of ions in 
the SHM is characterized by a large fraction of fast particles. 

In the case of a non-isothermic plasma (ZHM), the amplitudes of modes of the LF spec¬ 
trum (ion-sound waves) are of the same order in a wide range of wave numbers. In a cold 
plasma (SHM), the long-wave oscillations dominate in the LF spectrum. The energy of the 
LF held is found to be much lower than the total kinetic energy of ions for all the cases 
discussed above. Reducing the energy of the LF held with time happens due to the energy 
transfer to ions. 

The decrease of the damping rate of the HF held corresponds to the slowing of the 
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HF field burnout in the cavities and leads to the broadening of the HF spectrum that 
causes the deepening of the cavities and increase of the kinetic energy of ions ejected from 
them. It should be noted that as the absorption rate of the HF held decreases, the ion 
velocity distribution function approaches the Maxwellian distribution in both models under 
consideration. In a cold plasma, the intensity of the long-wave LF oscillations is high and it 
increases with the decrease of the absorption of HF modes. It is important to note that the 
total energy transferred to ions increases as the absorption of the HF spectrum reduces. 
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